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Abstract 

Recent observations of the TeV 7-ray spectra of the two closest active galactic nuclei 
(AGNs), Markarian 501 (Mrk 501) and Markarian 421 (Mrk 421), by the Whipple 
and HEGRA collaborations have stimulated efforts to estimate or limit the spectral 
energy density (SED) of extragalactic background light (EBL) which causes atten- 
uation of TeV photons via pair-production when they travel cosmological distances. 
In spite of the lack of any distinct cutoff-like feature in the spectra of Mrk 501 and 
Mrk 421 (in the interval 0.26 — 10 TeV) which could clearly indicate the presence of 
such a photon absorption mechanism, we demonstrate that strong EBL attenuation 
signal (survival probability of 10 TeV photon < 10~^) may still be present in the 
spectra of these AGNs. This attenuation could escape detection due to ambiguity of 
spectra interpretation between intrinsic properties of the sources and absorption by 
EBL. By estimating the minimal and maximal opacity of the universe to TeV 7-ray 
photons, we calculate the visibility range for current and future 7-ray observatories, 
and show that the Whipple 7-ray telescope should be able to detect (in 10 hours 
at a 5o" confidence level) a BL Lac object with properties similar to Mrk 501 dur- 
ing its peak activity located at distances up to z = 0.12. The proposed atmospheric 
Cherenkov telescope array VERITAS should be able to see such an object at least as 
far as z = 0.3. Finally, we show that the proposed experiments, VERITAS, HESS, 
and MAGIC, may even be able to actually measure the EBL SED because their 
observations extend to the critical 75 — 150 GeV regime. In this transition region 
a distinct "knee-like" feature should exist in the spectra of blazars, which is invari- 
ant with respect to their intrinsic properties. The change of the spectral index and 
flux amplitude across this knee, if observed for several blazars, will provide missing 
pieces of information needed to measure EBL in the wavelength range 0.1 — 30 /xm. 
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1 Introduction 



It was noted first by Gould and Sclireder (|T|) in the late sixties "that observa- 
tions of cosmic photons in the region lO'^^ to 10^^ eV would be of great value, 
since in this region absorption due to the cosmic optical photons is impor- 
tant. In fact, this may provide a means of determination of the optical photon 
density and of testing cosmo logical models." Although the authors pointed 
out that quasars may be the potential sources for such observations, and the 
ground-based Cherenkov technique could be the means of study, the real possi- 
bility to utilize this idea came one quarter of a century later with the discovery 
of TeV extragalactic 7-ray sources by the Whipple collaboration. These are 
low-redshift BL Lac objects Markarian 421 {z = 0.031) @, Markarian 501 
{z = 0.033) (D, and 1ES2344+514 {z = 0.044) (|). The detection of a more 
distant X-ray BL Lac object, PKS 2155-304 {z = 0.117), has been recently 
reported by the Durham group (|^). 



Interest in the absorption of extragalactic 7-radiation in recent years has been 
supported by two reasons: understanding the observational range for extra- 
galactic sources of the very high energy photons (above ~ 1 TeV) and estimat- 
ing the density of the EBL photons of rather low energy (below ~ 1 eV). The 
first motivation is of great importance for future 7-ray experiments, such as 
VERITAS (D, HESS (0), and MAGIC (g), and for testing of blazar emission 
models (|^; |10| 11; 12). The second attracts great theoretical attention, since 
EBL in the energy range ~ 10^^ — 10 eV contains a wealth of information 
on the early history of the universe providing integral constraints on galaxy 
and star formation, the initial distribution of stellar masses, metal and dust 
production, and the rate of re-processing of starlight to infrared wavelengths 



by the dust (|T3| ; |Tj). It is also possible that a significant fraction of the EBL 



is contributed by brown dwarfs, accreting black holes in active galactic nuclei 
(AGNs), and the decay of relic particles (|l^; |16D. Because the EBL contains 
a significant fraction of all energy released since the recombination epoch, its 
determination may set a limit on the energetic balance of the universe and 
constrain particular energy sources (ITBp. 



Diffuse Infrared Background Experiment (DIRBE) and Far Infrared Absolute 
Spectrometer {FIRAS) instrument on board the Cosmic Background Explorer 
{COBE) spacecraft (|l3) have been used to complete an effort to directly de- 
tect an isotropic EBL in ten photometric bands from 1.25 to 240 fim (|18|), 
and 125 — 2000 fim wavelength interval (plQl) . The detection of the EBL for 
wavelengths shorter than 140 /im have been greatly impaired by the pres- 
ence of strong backgrounds from interplanetary dust scattering and emission, 
stars and galactic foregrounds, and from interstellar dust emission. Removing 
these various contaminating radiations to detect a directionally isotropic sig- 
nal, which is presumably a unique feature of EBL, is a difficult task (EDI; pTI). 
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The COBE collaboration reported a detection of the EBL signal at 140 and 
240 /im, but established only upper limits in the region 1.25 — 100 /xm (18), 



the most important interval for attenuation of extragalactic TeV 7-rays. Thus, 
observation of TeV and sub- TeV extragalactic sources is currently the only rel- 
atively direct method which may provide new insight on the distribution of 
EBL in this wavelength range. 



There are two alternatives which can be used to constrain the EBL spectral 
energy density from 7-ray observations. First is to model EBL spectra based 
on the knowledge of galaxy and star evolution, and use it to find the opti- 
cal depth t{E,z) and survival probability exp{—T{E, z)) for 7-rays of energy 
E traveling from a source at redshift z. The spectrum of a particular source 
can be interpreted then as modified by energy- dependent EBL attenuation. 
Such predictions of the opacity of the universe to ~TeV photons have been 



reported in (p2| He B3; 24). Many more models have been suggested which 



derive the EBL SED based on observational constraints and theoretical as- 
sumptions (for a summary and references see (p!3D). These EBL predictions, 
which can be easily converted into a 7-ray optical depth, are roughly classified 
in two categories: backward and forward evolution models. The former models 
adopt current observations of spectra of galaxies as a function of luminosity 
and translate them backward in time using data on the redshift evolution of 
galaxy emissivity. An example of such empirically based calculations is given 
by Malkan & Stecker (pS]). The latter models tend to simulate the evolution 
of galactic systems forward in time assuming initial parameters which do not 
conflict with current observational data, e.g. (^; ^). A clear advantage of 
these models is the possibility to test various scenarios, and, if the EBL is 
observed to have certain features in a spectrum, it may be possible to deduce 
information about cosmology, galaxy and star formation, as well as absorption 
and radiation of the light by dust. Both backward and forward evolution mod- 
els do not guarantee, however, that all potential contributors to the EBL are 
accounted for and therefore more direct constraints on EBL spectral energy 
density are desirable. Such an alternative, which attempts to unfold the EBL 
SED from spectral measurements of extragalactic sources, has been explored 
recently by Stanev & Franceschini, and Biller et al. (|2^; In this paper 
we follow an analogous strategy to establish limits on EBL from the spectral 
measurements of two AGNs, Mrk 421 and Mrk 501, and define applicability 
conditions of such an approach. 



The interpretation of the EBL limits derived from TeV 7-ray measurements 
is difficult because the EBL absorption cannot be completely isolated from 
the properties of the source. Although, this problem may be overcome in the 
future when more TeV and sub- TeV sources are discovered and the physics of 
the 7-ray emission from blazars is better understood, at present, this feature 
of the method remains a major source of uncertainty. The credibility of these 
derivations of the EBL limits from TeV AGN 7-ray measurements is still not 
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accepted ([T8| ; |T^). The claim of the first measurement of the EBL by de Jager 
et al. (29) derived from an apparent "cutoff" in the spectra of Mrk 421 at 4.1 
TeV (|30|) has not been confirmed. More recent measurements of the spectra 



of Mrk 421 and Mrk 501 (|3l| ; p2t |33|) are of much better quahty but the inter- 
pretation is not unambiguous. The Whipple group ( p2D find curvature in the 
spectrum of Mrk 501 but conclude that the spectrum of Mrk 421 is consistent 
with a simple power law up to 10 TeV. The curvature in the spectrum of Mrk 
501 is confirmed by the HEGRA group with high precision (^3]). Because of 
the apparent differences in the curvature in the spectra of the two AGNs (sig- 
nificant only at the ~ 2 sigma level), Krennrich et al. ( pT]) suggest there is no 
evidence for EBL absorption. Stecker has recently revised an earlier predic- 
tion (^) and now suggests that the new spectrum of the flaring Mrk 501 shows 
EBL absorption and agrees with one of his EBL models (^If). The recent inves- 
tigation by Stanev & Franceschini of the Mrk 501 spectrum ( ^^ has indicated 
upper limits on the EBL which almost come into conflict with observational 
evaluations based on deep surveys of extragalactic sources. These authors con- 
clude that "if spectra at TeV energies for extragalactic 7-ray sources like this 
for Mrk 501 are confirmed with improved statistics, we may be forced to con- 
clude that the process of 7 — 7 interaction in the intergalactic space is more 
complex than expected and that the average intergalactic magnetic field is 
extremely weak {B < 10"^^ G)." 



In this paper we strive to reconcile these differing interpretations, we attempt 
to limit the scope of the potential constraints which 7-ray observations may set 
on the SED of EBL, and we try to confirm some of the earlier results derived 
from TeV 7-ray measurements. The structure of the paper is as following. The 
next section, §2, is a summary of the direct EBL experimental detections, up- 
per and low limits. These are important for unfolding SED of EBL from the 
spectra of blazars because they can be used to normalize the EBL contribution 
to 7-ray absorption and establish upper and low limits on TeV 7-ray survival 
probability. Then, in section §3, we consider details of the pair production 
process, which absorbs high energy photons, and point out several features 
important for understanding of this fundamental interaction. We also repeat 
why one would expect a cutoff in the spectra of TeV extragalactic sources 
caused by the EBL. The next section, §4, analyzes the recently reported spec- 
tra of two AGNs which provide data over a wide dynamic range, with high 
statistical quality. Although no apparent cutoffs have been detected up to 10 
TeV, both spectra are consistent with a non zero curvature. In section §5 we 
demonstrate that although a cutoff feature is possible, it is not necessarily the 
only indicator of EBL absorption and the lack of a cutoff in the spectra of 
AGNs does not automatically mean low density of EBL. Such an alternative 
is explored in detail in this section, and we derive a SED of EBL which eludes 
detection due to ambiguity in the interpretation of spectral measurements. 
Then we determine upper and low limits of this EBL component. In section 
§6, we use these data to predict maximal and minimal attenuation in the spec- 
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tra of a "Mrk 501-like" blazar at high redshift if such a hypothetical object 
demonstrates a peak of activity, and compare its flux with the sensitivities 
of the existing Whipple (^) and proposed VERITAS (|^) 7-ray observatories. 
The last section, §7, is devoted to a discussion. 



2 Summary of direct EBL detections, upper and low limits 



Fig. 1 depicts the current observational results on the EBL from UV to far- 
infrared wavelengths. The range from 1.25 /xm to 240 /im is covered by DIRBE 
results obtained after a complex background subtraction (|18]). Detections are 
found at 140 /^m and 240 /xm. Upper limits are reported for the other bands 
where foreground emission from the galaxy and interplanetary dust did not 
allow estimation of the EBL SED. A tentative detection has also been reported 
at 3.5 fim (|36D, derived with the use of the DIRBE upper limit and a lower 
limit from ground-based K-band galaxy counts (0). The latter, shown with 
the filled star in Fig. 1, was assumed to be a detection in this analysis. The 
upper limits in UV and optical are derived from all-sky-photometry measure- 
ments (|38| ; |39| ; ^T]) . For a more extended review of earlier results in this 



energy range see (^21). Dwek et al. (|T3|) also indicate recent detections in this 
wavelength range by Bernstein et al. (not shown). 

The lower limits obtained from galaxy counts have been derived from IRAS 
data (^31) for 25, 60, 100 /im. The two data points at 6.7 and 15 /im are deduced 
from ISO data ( ^4[) as reported in (0). Pozzetti et al. published ground- 



based measurements in K-band ( ^Sj) , and derived recently lower limits from the 
Hubble Deep Field (HDF) observations (^6]) covering the range 0.36 — 0.81 /im. 
The ground-based UV lower limit at 0.2/im is reported by Armand et al. (^T]). 
The solid line in Fig. 1 indicates a FIRAS detection of the extragalactic far- 
infrared background (|TP|) in the range 125 — 2000 fim. The cosmic microwave 
background dominates EBL for wavelengths above ~ 400 /im. 

The dotted line in Fig. 1 shows an example EBL prediction modeled by Pri- 
mack et al. ([1^) for a cold dark matter cosmology with non-zero cosmological 
constant and Salpeter stellar initial mass distribution (^Sj). The main features 
of the EBL curve, common for the majority of models, are a peak in SED distri- 
bution in the range from 0.1 to 10 fim formed mainly by starlight accumulated 
and redshifted through the history of the universe, the so-called stellar EBL 
component, and a second peak in the range 10 — 1000 /xm formed by dust re- 
emission of the previously absorbed starlight, the dust EBL component. The 
excess of SED prediction in the UV region can be adjusted to observations 
in this model either by a higher starlight extinction rate by the dust, or by 
changing the stellar initial mass distribution to the Scalo model (^91), which 
has less mass concentrated in high-mass stars and consequently produces less 
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UV light. The EBL predictions by the different models vary, and some of them 
generate less pronounced features, e.g. Malkan & Stecker calculated an EBL 
with a substantially smaller valley in the 10 /xm region (pSj). This difference is 
discussed in the last section. 



3 The absorption of extragalactic TeV photons 



We are interested in the process 7 + 7 
section is 



e"*" + e , for which the total cross 



(^(q) = o f{q) 



f{q) = <i 



(l + ,-|)lni±^-(l + ,)^ 



^ Ee l-cos(^)' 

where ctt = 6.67 x 10~^^ cm^ is Thomson scattering cross- sect ion, me = 0.51 
MeV is the mass of the electron, E and e are the high and low energies of 
the photons, and 6 is the collision angle (see, for example, (|50|) ). The /(g) 
function, shown in Fig. 2, reaches its maximum at g = 0.508, indicating that 
for a head on collision the peak of the interaction cross-section of the 7-ray 
photon of energy ~ 1 TeV corresponds to pair production with a soft photon 
of energy e ~ 0.5 eV (A ~ 2.5 /im) 



A = 2.41/im 



E 
TW' 



(2) 



The absorption probability of the high energy photon per unit path length d/ 
by isotropic diffuse radiation with spectral density dn{e)/ de, is given by (see, 
for example, original calculations by Gould & Schreder (|5T1)) 



dr 3 /",/■, ,,^1 Ix 



Ee 



1 — COS 9 



X - 



dn{e) 
de 



(3) 



Eq. (3) can be transformed to the following form 
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dr 3 7 , dn(e) ^,ml^ 

ml 
1 

F{q) = 2q^ J f{x)x-^dx. (5) 

Therefore, attenuation of the high energy 7-rays by isotropic background pho- 
tons peaks at the maximum of F{q) {q = 0.28), shown in Fig. 2. This gives 
a characteristic energy of the soft photon £ ~ 0.9 eV for the most effective 
attenuation of a 1 TeV 7-ray 



A = 1.33//m;^. (6) 

This energy is almost a factor of two higher than for a head on coUision given 
in Eq. (2). The upper axis in Fig. 1 indicates energies of the 7-ray photons 
coupled to the wavelengths of the EBL photons via Eq. (6). 

The lowest energy of a soft photon important for interaction with a hard 
photon of energy E is determined by the threshold for pair production in a 
head on collision 



A = 4.75//m;J^. (7) 

Thus, for example, a 100 TeV 7-ray begins to feel the presence of the CMB, 
but it is attenuated mostly by the dust EBL component. 

The highest energy e relevant to interactions with extragalactic photons of 
energy E should be derived from the properties of the function F(q). The 
asymptotic behavior of this function at small q (large e) 



lim F( — ) 2 — (In — 2) 

£-»oo ^Es' Ee^ ml ' 

converges very slowly to zero, such that the integral 



2 

E 

is divergent at the upper limit. This indicates that interaction of a hard photon 
of energy E extends far below the characteristic wavelength given by Eq. (6), 
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and suppression of the integral (4) divergence at large e is due to the behavior 
of the EBL spectral density (\n{e)/ de. It follows then, that the SED of EBL 
itself defines the characteristic energy of the soft photon and range of the 
photon's energy important for intergalactic absorption of a 7-ray with energy 
E. 

As an important example we consider a SED of the form 



dn{e) a f e \ ^ 



de eV VeV 
and derive the differential absorption probability 



2.Ta3(/3)(4f" , (10) 

\ TJl/fi J J 



Depending on the power index /?, the main contribution to this integral comes 
from the different regions g, and, therefore, from different energy ranges of the 
EBL field. For (3 = 2.55 the function q^'^^F{q) is shown in Fig. 2. It peaks 
at g = 0.41 and 90% of the contribution to the integral g{2.55) = 1.47 is 
accumulated when q changes from 0.13 to 0.80. Only 5 percent tails are left 
below and above this interval. Thus, the range of soft photon wavelengths 
important for attenuation of 7-ray of energy E is 



dr _ 
'dl~ 

9m= 



0.62/im < A < 3.8/im (11) 
The similar limits derived for /3 = 2.00, when value of ^'(2.00) = 1.25, are 



0.37/im < A < 3.6/im 

It can be seen that while low energy 
the same for these two examples, the 
1.7 extending the range important for 
energies in the case P = 2.00. 



(12) 

end of soft photons remains almost 
high energy end changes by a factor 
attenuation of 7-ray photons to high 



For a 7-ray of energy, E, which travels cosmological distances from a source 
at redshift z we generalize Eq. (4) following, for instance (^2]), to obtain the 
optical depth 
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t{E,z) = ^ [ VTT^ dz J de (13) 





no ho 8 Ho' 



where Ho = 100 ho km Mpc~^ is the Hubble constant, 0.5 < ho < 0.85 is 
the normahzed Hubble expansion rate, c is the speed of light, no = 4.33 x 10""^ 
cm~^ is the characteristic density, and dn{e)/ de is the present-day spectral 
EBL density. It is also assumed that the density parameter of the universe 
Q = 1, cosmological constant Q/^ = 0, and evolutionary considerations of the 
EBL spectral density are neglected, which is justified for {z < 0.3), implying 
that the EBL has been produced mostly at higher redshifts (p^). 

The appearance of a cutoff in the spectra of TeV extragalactic sources has 
been suggested in one of the pioneering works on TeV 7-ray absorption by 
Stecker, de Jager and Salamon ( |5^ ) for quasar 3C 279 {z = 0.54) and later 



several predictions (^; |2^; |56|; ^T]) for Mrk 421 have been published. If the 
differential energy spectrum of a source obeys a power-law, characterized by 
index 6 and flux J one would expect to observe a spectrum modified by EBL 
absorption 

^ = J E-'eM-r{E,z)). (14) 
For an EBL spectral density (9) and small z, the optical depth is given by 



\a g{(3} zj 

and the prediction for an observed spectrum of the source is then 



If (3 is large (as in the original paper which claimed possible EBL measure- 
ment ( PI?D from TeV 7-ray observations: (3 = 2.55, 6 = 2.1), a sharp cutoff must 
exist in the spectrum because the power-law term in the exponent dominates 
the logarithmic term at energies above E^.. A lack of events in the preliminary 
spectrum of Mrk 421 in bins above 3.4 TeV and 5.1 TeV was interpreted as 
such a cutoff providing an estimate of ii^c = 4.1 TeV. Due to a rapid relative 
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change of the two terms in the exponent, the logarithmic term dominates for 
energies below 0.73 = 3.0 TeV leaving almost no trace of the cutoff at 
lower energies (< 0.3 Ec). Although this detection has not been confirmed in 
later observations (see next sections §4) the approach to place limits on EBL 
by "non-detection" of the cutoff feature in AGN spectra remains legitimate. If 
the spectral measurements of the source extend without a cutoff to an energy 
Ec, the EBL spectral density (9) is constrained by 



ho /TeV\ 



(3-1 



g{(3) z\ Ec J 



(17) 



It is important to note that as E^. changes in observations, both the a limit 
and the "window" of EBL photons, which did not produce a cutoff feature at 
Ec, change according to Eqs. (17, 11, 12). Because this limit depends on the 
EBL SED model, the model itself must be justified by current theoretical con- 
siderations in the corresponding wavelength interval. For example, the region 
10 — 100 /im is believed to be approximated well with a SED of the form (9) 
with (3 in the interval 2 — 3. Such approximation, however, is invalid for the 
1 — 10 /im wavelength range. Due to the sensitivity of the 7-ray absorption 
rate to a rather extended interval of EBL photons (properly explained in Eq. 
(8)), the optical depth t{E, z) is not particularly sensitive to small details of 
the EBL distribution, they are smoothed out by the integration (Eq. 13). On 
the other hand, the absorption effect on 7-rays may be exponential, at least 
for some EBL models (Eq. 16); this is the feature which allows one to establish 
EBL upper limits by a non-detection of the cutoff. Recent Mrk 501 observa- 
tions by the HEGRA collaboration ( P^D seem to indicate a rapid decline in 
the spectrum of this AGN above 10 TeV. If the systematics of HEGRA ob- 
servations in this energy regime are well understood and these measurements 
are confirmed, an interpretation of this decline in the Mrk 501 spectrum by 
EBL absorption may constrain the spectral density in the 10 — 70 yum range. 



4 Mrk 421 and Mrk 501 spectra 



Krennrich et al. (|32D have recently published the compiled TeV spectral data 
of two AGNs, Mrk 421 {z = 0.031) and Mrk 501 {z = 0.033), shown in Fig. 3. 
The high statistical quality of these measurements and wide energy range, 
covering 0.26 — 10 TeV, allowed for the first time detection of curvature in 
the Mrk 501 spectrum during its high state of activity. This variable blazar 
shows changes in the TeV 7-ray flux on a time-scale of several hours ([58|) 
with no evidence for temporal variability of the spectral shape observed by 
the Whipple ( |59D and HEGRA (^) collaborations. The spectrum remains 
statistically invariant from the absolute flux which may increase up to several 
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Crab units0 during flaring activity ( pSf |60|). 

Mrk 421 is also highly variable showing bursts of radiation on a sub-hour time 
scale (^2j), with flux variations of up to 10 Crab units (^). The spectrum of 
this blazar was derived from three data sets with average fluxes of 2.8, 3.3, and 
7.4 Crab units. Neither of these sets has revealed any statistical differences in 



the spectrum shape justifying compilation of the data (|6^) 



Both spectra are rather featureless in the energy range 0.26 — 10 TeV and do 
not show an apparent cutoff which could be interpreted as an EBL absorption. 
Both spectra are well described by a flux dependence 



= J exp (-51n(E/TeV) - c\n\E/TeV)) , (18) 

where 5 is a spectral power index, and c is a curvature. The best flt of such a 
representation of the data is shown in Fig. 3. It can be seen from the flgure 
that the spectra differ; a test places the chance probability that they arise 
from the same parent distribution as 4 x 10"'^ (^2]). The differences in the two 
spectra must be due to differences in the intrinsic emission spectra since they 
are at practically the same redshift; however the fraction of curvature that is 
common to both spectra could be explained by EBL absorption. 

The differences can arise from differences in the spectral index, 6; in the cur- 
vature parameter, c, or in a combination of both. The simplest hypothesis is 
that a spectrum can be flt by a simple power law (c = 0) and, as discussed 
in (0), the spectrum of Mrk 421 is compatible with this hypothesis. The 
spectrum of Mrk 501 is not. The spectral indices of the two blazars are sta- 
tistically different. Fig. 4 depicts the chance probability that a certain value 
of spectral index can be associated with the data on each object. The curva- 
ture parameter in these estimates has been flt to the best value at a given 
index. The statistical quality of the data is enough to rule out the hypothesis 
of equal spectral indices for the two AGN. A similar analysis of the curvature 
parameter is not deflnitive. Although the best flts to c for each source are 
different (Mrk 421 0.12 ± 0.04; Mrk 501 0.20 ± 0.03 (0)) it is also possible 
to flt both spectra with the same curvature parameter. Fig. 5. Some or all of 
this curvature could arise from EBL absorption. The peak of the probability 
{p^2 = 0.1) is at c = 0.17 with the chance probability less than 10~^ that the 
curvature parameter exceeds 0.3. 



In the Whipple paper (|32D the differences in the spectra of two AGNs were 
interpreted in the context of synchrotron-inverse Compton (SIC) blazar emis- 
sion model, and it was concluded that the harder spectral index of Mrk 501 



^ The Crab unit is defined as the 7-ray flux above 300 GeV from the Crab Neb- 
ula (|l|), i.e. N^iE > 300 GeV) w lO'^ m-^g-i. 
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and its non-zero curvature is not inconsistent with such a model, if X-ray 
observations of these objects are attributed to synchrotron emission. It was 
pointed out that the difference in the energy spectra of Mrk 421 and Mrk 501 
must be intrinsic to the sources and not due to intergalactic EBL absorption 
since both AGNs are located at the same distance z = 0.03. It is shown in 
Fig. 5 that a substantial part of the curvature is not ruled out as common for 
both blazars, and therefore difference in their curvatures is to be considered 
as tentative only. Moreover, explaining the curvature of Mrk 421 and Mrk 501 
as dominated by EBL absorption seems plausible in the context of SIC model 
because it is easier to understand then the observed spectral shape indepen- 
dence of these blazars relative to their absolute flux during episodes of activity. 
When the intrinsic spectrum is almost a pure power-law, it remains invariant 
with respect to both amplitude and energy re-normalization. If a substantial 
curvature is present in the intrinsic spectrum, the latter invariance is lost, 
severely constraining models of blazar flaring (see for example (|66D). If, how- 
ever, the curvature of Mrk 501 is mostly due to EBL absorption, which is 
not inconsistent with observations of Mrk 421, then changes in the intrinsic 
hardness ratio of the AGN during flaring may be partially masked by this 
effect. 

In the recent submission (^) by the CAT collaboration on the spectral prop- 
erties of Mrk 501 the Whipple result on spectral index and curvature of this 
object has also been confirmed (2.24, 0.22). Unlike the Whipple and HEGRA 
observations, CAT data seem to indicate spectral variability depending on the 
state of emission. The detected change of hardness ratio is interpreted in the 
frame of a SIC blazar emission model with synchrotron target photons, and 
an interpretation of Mrk 501 flaring activity is suggested as being due to injec- 
tion of a population of energetic electrons and/or an increase in acceleration 
efficiency. Although the non-detection of spectral changes by the HEGRA and 
Whipple collaborations is explained by high energy threshold and low statistics 
respectively, we think that this result needs further confirmation to reliably 
disentangle the dynamic and static parts of the spectral index and curvature 
parameters. The latter should presumably provide a better constraint on the 
EBL absorption effect. 

Both the HEGRA and CAT collaborations have detected a signal above 10 
TeV in the spectrum of Mrk 501 during its highstate activity (|6^; in 1997. 
The spectral measurements performed with the HEGRA stereoscopic system 
indicated the presence of an exponential cutoff at 6.2 TeV in the power-law 
spectra of this object. In the region 0.5 — 10 TeV, however, the reported 
spectrum is consistent with the Whipple result and can be well described 
by the power-law with a curvature term. The region above 10 TeV does not 
follow this dependence and shows a substantially more rapid decline in the 
AGN flux. Although the cutoff feature on its own does not allow a unique 
interpretation in terms of EBL limits and/or physical processes generating 
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radiation at the source, it provides nevertheless a set of interesting constraints 
suggested in (68), among which an increase in the SED of EBL in the region 
30 — 70 /im is a possibihty. 



Because the cutoff above 10 TeV in the spectra of AGNs has been observed 
only for one object and only by one instrument we do not discuss this fea- 
ture in our paper due to a variety of alternatives in its explanation. Instead, 
we refer the reader to several recent publications (|69| ; ffOt ffll) . Alternatively, 
the region of the Mrk 501 spectrum below 10 TeV is well established and 
confirmed by several groups including Whipple, CAT, HEGRA, the Telescope 
Array (0), and the Crimean observatory (|73D . In addition, spectral measure- 
ments of another AGN at a similar redshift, Mrk 421, are known in the same 
energy range and both spectra a well fitted by a power-law with the curvature 
term. The HEGRA collaboration has reported recently a Mrk 421 spectrum in 
a low state of its activity (^) during 1997-1998. Likely due to limited photon 
statistics, no nightly variability has been found in the spectra of this blazar. 
The time averaged energy spectrum has been found to be very steep with an 
index of 3.09 ± 0.07 in the energy range 0.5 — 7 TeV. Although the signifi- 
cance of the difference between Whipple and HEGRA measurements and its 
impact on flaring spectral variability of this source remain to be understood, 
in both measurements the time averaged AGN spectrum is consistent with a 
pure power-law, but an exponential cutoff or curvature term can be equally 
well fitted to the data. We use Whipple data which spans the same energy 
interval (0.26 — 10 TeV) to conservatively conclude that curvature, c, of the 
spectra of the two AGNs caused by intergalactic absorption does not exceed 
0.3. The curvature of both AGNs may be due to a common cause, if so, its 
most probable value is ~ 0.17. 



5 EBL unfolding 



The detailed investigation of the spectra of two AGNs has not revealed any 
"sharp" feature. The non-detection of the cutoff, however, does not automat- 
ically mean low EBL spectral density and negligible intergalactic absorption 
of TeV 7-rays. The featureless spectra of blazars, in fact, allows a non-trivial 
interpretation. The underlying reason for this is that EBL absorption is diffi- 
cult to separate from the properties of the source. It can be seen from Eq. (16) 
that if parameter (3 of the EBL SED is close to 1, the intergalactic absorption 
does not produce a distinct feature in a source spectrum. We investigate such 
a possibility in this section. Our strategy will be just opposite to a prediction 
of optical depth from EBL spectral density. Instead we accept that r(i?, z) 
can be approximated by 
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r(i?, = Co{z) + c,{z) In (;^) + c,{z) In' (^) , (19) 

which is supported by observations in the energy interval (0.26 — 10) TeV for a 
small z. One can think of such an approximation as an expansion series of the 
ln(£^/TeV) variable, when higher orders produce undetectable contributions in 
a given energy interval for known spectra of blazars. Using this optical depth 
representation, we unfold the corresponding EBL spectral density. 

To simplify our considerations, we assume at the moment that the redshift of 
the source is very small {z <^ 1). Using formula (13) we obtain an integral 
equation 



J^s^nfy-^riE, (20) 

E 

which can be resolved with the use of a Mellin transformation. Omitting this 
derivation, we give the final result 



e 



dn{6) riQ ho f j^,Tn1, , , dE 



de z 

™2 



cr+ico 
a—ico 

1 

G{s) = J q'-' F{q) dq. 



It can be seen from Eq. (21) that the formal solution dn{e)/ de which satisfies 
the optical depth (19) may be represented as 

iM£)=!^2L„_p,l,± + p,l„a±) (22, 
de e \ £o £o/ 

with characteristic energy Eq and some constants Po,Pi,P2- We use this form 
of the EBL spectral density to connect (po,Pi,P2) to (co,Ci,C2) of the optical 
depth, conducting a direct calculation of integral (13) with arbitrary, but not 
large redshift {z < 0.3) required by an applicability condition of this formula. 
One can find that 



1 ^ 

Ci^j-Y.Pj'^jii^)^ (23) 

ho ,=1 
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where {i = 1, 2, 3), and tji{z) is a redshift dependent matrix. The introduction 
of a finite redshift, z, does not change our conclusion about the form of the 
EBL SED (Eq. 22) since the integration order in integral (13) can be changed, 
and z should be considered then as a parameter. We also use a characteristic 
energy sq = 9.3 eV (A = 0.13 fim) of the isotropic EBL which corresponds to 
a peak interaction with 100 GeV 7-rays and truncate the SED at this energy. 
The current TeV and sub-TeV measurements are not very sensitive to the 
EBL distribution in UV region. The truncation, however, will introduce an 
error in the optical depth for low energy 7-rays [E ~ 100 GeV) and larger 
z, limiting the applicability of our considerations to the 0.1 — 10 TeV energy 
regime. We expect, though, that attenuation of 100 GeV 7-rays will be rather 
small for the range of redshifts investigated here. 

Three examples of EBL spectral density described by Eq. (22) are shown in 
Fig. 7. It is a curious coincidence that these three parametric SED of the 
EBL can describe well the expected peak in the 0.1 — 10 fim wavelength 
range of the EBL energy density produced by starlight emitted and redshifted 
during evolution of the universe. There is no a priori reason why the current 
observational window of TeV astronomy (0.26 — 10 TeV) and requirement 
that EBL absorption of 7-rays must produce almost power-law modifications 
to the spectra of extragalactic sources would necessarily require that the EBL 
SED, which is allowed to have a peak in the wavelength range 0.1 — 10 /im, 
would coincide with the theoretical expectation based on completely different 
assumptions. In other words, if the starlight peak were to be shifted to a longer 
wavelength region by a factor of 5 — 10, any reasonable amount of EBL would 
necessarily produce a "knee-like" feature in the (0.26 — 10 TeV) spectra of 
the sub- TeV sources. The starlight EBL peak is located in such a wavelength 
region that it eludes its apparent detection in the TeV and sub- TeV spectra of 
blazars due to ambiguity in the interpretation of the spectral measurements. 

The direct upper and low limits on the EBL SED (see Fig. 1) restrict a possible 
choice of parameters po, pi, and p2- For example, pq, which defines the SED 
of EBL at A = 0.13 ^m, is likely limited to (0.025 — 0.45) range by Armand 
et al. ( P7| ) 0.2 /im lower and by Bowyer (^) 0.165 /im upper limits. While 
the exact choice of this parameter affects the optical depth of 7-rays with 
energies around 100 GeV at large z, it has a rather negligible effect on the 
attenuation of photons in the 0.3 — 10 TeV energy range. Thus its value is not 
constrained by current TeV observations of the low redshift sources. With the 
reservation that we do not predict correctly a small optical depth of ~ 100 
GeV photons we fix this parameter somewhat arbitrarily, at 0.025, assuming 
that EBL attenuation will unlikely be lower than the one we find. 

The limitations imposed on the choice of parameters pi and p2 are shown in 
Fig. 6. The most severe constraints are due to the DIRBE upper limit at 3.5 /im 
and due to lower limits found by Pozzetti et al. (^) for K-band and by Oliver 
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et al. (1^) at 15 /im. The upper limits in the wavelength range 0.165 — 0.512 
/im may potentially constrain the pi parameter as it is shown. At the same 
time these limits affect little the value of p2 which mostly determines the EBL 
SED at longer wavelengths and absorption of TeV 7-rays. The pi parameter 
may affect the optical depth in the range 0.3 — 1 TeV. However, if survival 
probability of such photons is close to one, formula (22) may not be well 
justified and in this case limits on pi from the EBL UV region may become 
invalid. On the other hand, if the opacity of intergalactic media to such 7-rays 
is large (> 1), these limits must be taken into account. Gamma-ray spectral 
measurements of sources with larger z, such as the recently discovered (^) 
X-ray BL Lac PKS 2155-304 {z = 0.117), may provide important information 
to resolve this question. 

In the case of Mrk 421 and Mrk 501, the matrix tji{z) which connects the 
parameters of the EBL SED to parameters of optical depth is 



^0.043 0.007 -0.002^ 



tji(0.032) 



0.080 0.041 0.003 
0.182 0.168 0.041 



(24) 



The elements of the tji upper triangle are not exactly equal to zero because 
we truncated the SED of EBL at 9.3 eV, but they are negligible. One can see 
that curvature in the spectra of these AGNs is driven by the p2 parameter 



C2 = — 0.041p2. (25) 
ho 

Although the upper limit on C2 < 0.3 (see §4) is consistent with the DIRBE 3.5 
/im limit, it provides no additional constraint on p2 < 7.3/io for a reasonable 
range of the normalized Hubble constant (0.5 < ho < 0.85). The preferred 
value of C2 = 0.17, if considered to be dominated by EBL absorption, may 
give an estimate of p2 = 4.15/io which is not ruled out for any ho from the 
indicated range. 

The change of the AGN spectral index due to EBL absorption 



ci = — (0.041pi + O.I68P2) (26) 
ho 

can range from 0.12//io to 0.61//io- This is a large interval since it is possible 
to have (~ 0.2 — 1.0) EBL contribution to the 2.28 spectral index of Mrk 501. 
Currently we are not able to use this EBL constraint. But once the mechanism 
of blazar emission is understood and the Ci contribution is isolated or limited 
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from the measurements of the AGN spectrum, then Eq. (26) may become an 
important restriction on the EBL SED. 

The total flux from a blazar is another parameter subjected to ambiguous 
interpretation. We find that the optical depth for a 10 TeV photon is given by 

r(10 TeV, z = 0.032) = -^(O.OSpo + 0.19pi + O.TQpa)- (27) 

ho 

The line of equal survival probability of such a photon is shown in Fig. 6. The 
maximal and the minimal optical depth allowed in the limited region of (pi, P2) 
parameters is reached at (0.,3.6) restricted by the DIRBE 3.5 fim point, and 
almost at (0., 0.7) constrained by the Pozzetti et al. ( P5| ) 2.2 /im measurement. 
Thus the possible range of optical depths is given by 0.55 /ho — 2.88//io or 
approximately from 1 to 5. In the latter scenario survival probability for such 
a photon is less than 1%. Fig. 7 shows the EBL SED for these limiting cases of 
maximal and minimal attenuation of 10 TeV 7-rays. This figure also depicts 
an example of the SED of EBL which satisfies a tentative EBL detection point 
reported recently by Dwek & Arendt (|36|). 

Concluding this section we point out an important limitation which concerns 
our SED solution for EBL given by Eq. (22). This limitation arises from the 
general problem of "ill-defined" unfolding tasks being very sensitive to small 
fluctuations in the source function (0), in our case optical depth t{E). In 
principle, the EBL unfolding procedure given by Eq.(21) can be performed 
for arbitrary t{E). We, nevertheless, chose a smooth approximation (19), and 
obtained a smooth solution (22). If t{E) is slightly perturbed, even with sta- 
tistically small fluctuations, the unfolded SED of EBL may acquire large dis- 
tortions. This fact, mainly originating from the asymptotic properties of F{q) 
(see Eq. 8), simply indicates a deficiency in the complete reconstruction of 
the soft photon density by unfolding of the optical depth. We assume, there- 
fore, that the suggested EBL solution (22) describes a smooth component of 
the SED. Small perturbations or spectral line features with small amplitudes 
are possible in the EBL SED as long as they are not detectable in t{E) be- 
cause integration (20) smoothes them out drastically. We also note that if the 
EBL SED is very small, close to the indicated lowest limit, then formula (22) 
may also become inapplicable, except in the high energy end, because EBL 
absorption will be virtually undetectable in the rest of the spectrum. 



6 Attenuation of 7-rays by EBL 

An uncertainty in the interpretation of AGN spectral measurements when 
determining an attenuation of the AGN absolute flux and change of its spectral 
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index because of the EBL intergalactic absorption creates difficulty in the 
prediction of the visibihty range of TeV blazars. We use our maximal and 
minimal estimates of the EBL SED, found in the previous section, to predict 
the maximum and the minimum effect of EBL on the observable spectrum 
of a "Mrk 501-like" blazar at higher redshifts if such a hypothetical object 
demonstrates a peak of activity. Then we compare its flux with the sensitivities 
of the existing Whipple (^) and proposed VERITAS (^) 7-ray observatories. 

We extrapolate the optical depth r{E, z) to higher redshifts by calculating 
the tji{z) matrix. Fig. 8 demonstrates the dependence of co(-2), Ci(z), 02(2;) of 
the parameterized t{E,z) expression (Eq. 19). This figure shows three sets 
of curves marked 1, 2, and 3. The ffist set corresponds to a maximal allowed 
absorption. The third describes the minimal EBL effect based on the already 
resolved EBL SED from galaxy counts in deep surveys of extragalactic sources. 
The second set of curves indicates the SED of EBL defined by (^1,^2,^3) = 
(0.025, 0., 4.15/10); which is preferable if curvature in the spectra of Mrk 501 
and Mrk 421 is dominated by EBL absorption (see §4). 

Figure 9 summarizes changes in the spectrum of a "Mrk 501-like" blazar placed 
at redshifts z = 0.03,0.06, ...,0.3. The total flux of this hypothetical AGN is 



normalized to 10 Crab units, the highest flux ever observed from Mrk 421 (|62D. 
The sensitivities of Whipple and VERITAS, marked "W" and "V", are indi- 
cated (|^). They are derived for certain cuts on parameters of the atmospheric 
shower images and correspond to a 5a detection in 10 hours, or detection of 
at least five photons in the regime of low background. The sensitivity shown 
reflects the standard, small zenith angle, observation mode. In the high en- 
ergy end it can be substantially improved by utilizing the large zenith angle 
observation technique ([3^). Two plots "a" and "b" demonstrate the absorp- 
tion effect for the highest allowable and the lowest resolved EBL. The change 
of the flux at 100 GeV when z increases is mostly due to a simple distance 
factor ~ z'"^. The intergalactic absorption of such 7-rays is small relative to 
this effect, and our EBL attenuation prediction, which is not reliable in this 
energy regime, is masked by this factor. For the lowest EBL SED, the dis- 
tance effect dominates absorption for all shown curves in the energy interval 
below 2 TeV. In this case, the Whipple telescope should be able to detect such 
a blazar up to 2; ~ 0.18 during a few days long episode of activity. For the 
worst-case scenario of intergalactic absorption the Whipple visibility range to 
such a hypothetical object is limited to z ~ 0.12. Thus, it seems that the high 
variability of blazars is the reason which so far prevents us from detection of 
such more distant objects, or luminosity of potential TeV blazars is too low 
to identify them at high redshifts as X-ray emitters. The recently detected X- 
ray selected BL Lac PKS 2155-304 has been observed to have maximal X-ray 
output in the range 0.05 — 0.5 keV during active states (|75D, which is similar 
to Mrk 421 peaking at a few keV (ffGl). This object showed a similar observed 
luminosity in X-ray but it is located at z = 0.117. The detection of such a 
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more luminous but more distant BL Lac is certainly not limited to 2; < 0.12 
and may extend to higher redshifts. 

The new generation of 7-ray observatories, such as VERITAS (D and HESS (|^, 
should be able to detect "Mrk 501-like" blazars at least up to distances of 
z = 0.3 as shown in Fig. 9a. Due to their much improved ability to survey the 
sky, the list of identified TeV and sub-TeV emitters will substantially increase. 
It is interesting to denote then, that these proposed projects may indeed mea- 
sure the EBL because of the extension of their energy range to below 100 GeV. 
The SED of the EBL is expected to drop very rapidly in the UV above 2 eV 
making the universe mostly transparent to 7-rays with sub-100 GeV energy. 
Therefore, in the vicinity of 100 GeV one should see a "knee-like" transition 
region in the spectra of AGNs where one power-law is replaced with one with 
a larger spectral index. The position of this feature should be independent of 
characteristics of blazars, and the accumulated change of spectral index and 
change of flux amplitude for several objects at known redshifts will define the 
EBL SED as shown in our calculations. The more accurate EBL unfolding 
will be possible because measurements in the region 100 GeV and below pro- 
vide a reference point in the spectra where no EBL absorption effect should 
be present. The remaining uncertainty related to a particular source intrinsic 
spectrum, perhaps, will be eliminated when more sources are detected. The 
visibility range of VERITAS and HESS will not be limited to z < 0.3 if TeV 
and sub-TeV emitters can have a much higher luminosity than Mrk 501. 



7 Discussion 

The DIRBE collaboration has completed analysis of the data on direct mea- 
surements of the EBL spectral density and published detections and upper 
limits (|IBp, which in many wavelength bands substantially exceed limits in- 
teresting for the theoretical understanding of the history of the universe. Cur- 
rently, only TeV 7-ray astronomy can potentially produce any additional con- 
straints in this field. At the moment though, its role as a contributor to our 
knowledge about EBL was characterized by the DIRBE collaboration in a fol- 
lowing conclusion (0). "Currently, there is no evidence for any intergalactic 
absorption in the spectrum of the three TeV sources detected to date. Claimed 
detections of, or upper limits on the EBL from TeV 7-ray observations should 
therefore be regarded as preliminary at present." We acknowledge the fact 
that some papers seeking to use TeV measurements to detect or limit EBL 
in some cases may have been based on preliminary judgments of unverified 
results. However, reliable conclusions derived from TeV 7-ray observations do 
exist. 

The "no cutoff" approach to derive upper limits on EBL discussed in §3 is a 
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valid method as long as assumptions of the model are justified in the relevant 
wavelength interval. It might be arguable if the original calculations, presented 
in the papers ( |5^ ; p9|), have always satisfied this requirement. According to 
the current paradigm, though, the wavelength regime 10 — 40 /im is expected 
to be dominated by emission from dust, for which EBL spectral density (9) 
with P = 2.0 — 3.0 provides an acceptable description. The absence of a cutoff 
in the spectrum of Mrk 421 at least up to = 8.8 TeV (/? = 2.55), or = 7.0 
TeV (/5 = 2.00) translates to the upper limits on the SED of EBL (see Eqs. 
(17 11, 12) and Fig. 3) 



Pa < 1-6 X IQ-^ho eV cm'^ , 2.0 /im < A < 37 pm (28) 

/ A \ 

Pa < 1-0 X IQ-'^ho — eV cm'^ , 4.1 pm < A < 40 pm. (29) 



10 pm 



Unless the intrinsic spectra of both Mrk 421 and Mrk 501 are exponentially 
rising to exactly compensate EBL absorption, these limits constrain rather 
strongly existing EBL models in this wavelength range (see Fig. 10) as they are 
almost two orders of magnitude lower than the DIRBE upper limits at 12 and 
25 pm. Kashlinsky et al. ([77| ; [78|) derived comparable limits from the search for 
a fluctuation signal in the DIRBE maps which reflects the spatial correlation 
properties of the galaxies implying that EBL comes from the redshifted light 
emitted by them. This method relies upon an extended set of assumptions 
on the evolution of such parameters as the galaxy SED, galaxy luminosity 
function with redshift, galaxy clustering, dust absorption and re-emission. 
The effects of intervening dust, among these, are less understood, especially 
for this wavelength range because the EBL re-processing by dust could have 
occurred outside of the galaxy which produced it. The limits derived from 
the "no cutoff" observations of 7-rays, on the other hand, require only the 
assumption of the EBL SED behavior in a certain wavelength interval and 
they are applicable only to models which satisfy this condition. If one considers 
the other limits for this wavelength range (shown in Fig. 10) derived from the 
different assumptions explained below, one can conservatively limit the SED 
in 10 — 40 pm interval to be below 10 nW m~^ sr~^. 



No apparent cutoff due to intergalactic absorption has been observed in the 
0.26 — 10 TeV spectra of the two most-studied AGNs. However, the recently 
reported spectra of Mrk 421 and Mrk 501 ( |32D are not inconsistent statis- 
tically with a common major contributor to the curvature of the spectra of 
these blazars, perhaps EBL absorption. This explanation might be preferable 
because almost pure intrinsic power-law of these sources provides a larger de- 
gree of invariance to support the lack of spectral shape variability observed 
during flaring activity in Mrk 501. This result, however, needs further conflr- 
mation. 
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The fact that no distinct cutoff feature has been observed in the spectra of two 
blazars allows a non-trivial interpretation of the EBL SED, as shown in this 
paper. The spectral energy density should have a certain form and coincide 
with the expected peak in this wavelength region of the SED from starlight 
emitted and redshifted through the evolution of the universe. It is rather a 
peculiar coincidence that TeV 7-ray observations currently span an energy 
range in which the contribution from the EBL absorption to the spectra of 
extragalactic sources can be masked easily by unknown properties of the source 
which has a spectrum close to a pure power-law. Extending the energy range 
to either higher energies or lower ones should identify a cutoff and knee-type 
features generated by the dust and stellar components of EBL, respectively. 

At the moment the upper limit on the EBL SED can be found (see Fig. 10) 
based on two assumptions. First, we normalize the EBL spectral density at 
the 3.5 /im DIRBE upper limit and, second, we propagate the SED solution to 
shorter and longer wavelengths requiring that such an EBL SED may generate 
changes only in the source flux amplitude, spectral index of the power-law, 
and spectral curvature. These are allowed by observations of Mrk 501 and Mrk 
421. The expected rapid decline of the EBL SED in the UV region is also used 
in these considerations. Our limits in the range 1 — 10 /im are consistent with 
those previously derived by Biller et al. ( ]28|) and by Dwek & Slavin (|56D. The 
former work also used the DIRBE upper limit for normalization and propa- 
gated the allowed SED solution by limiting the range of spectral index changes. 
Dwek & Slavin were first to study 7-ray absorption by the stellar component of 
EBL and found high model dependent upper limits in this wavelength range. 
This is the first work which indirectly proved that a substantial amount of 
the EBL produced by redshifted starlight may escape detection in TeV 7-ray 
observations of a source with unknown absolute fiux and spectral index. 

The limits derived from TeV 7-ray measurements in 1 — 10 /im region seem 
to be in good agreement with each other except the very low ones given by 
Stanev & Franceschini (|2^). These authors also attempted to unfold EBL 
spectra indirectly from the measurements of Mrk 501, but the absolute fiux 
normalization and power-law spectral index were allowed to vary during the 
fitting procedure. As we showed in §5, an EBL SED component of the form 
e{j)Q — pi\\i{£ / eo)) with arbitrary constants is likely to be undetectable with 
such a method. Perhaps existence of this hidden invariance was the reason 
for unstable behavior of the solution which has been found to be crucially 
dependent on the manner in which power-law fiuxes are normalized. If our 
supposition is valid such a component can be added to the very low upper 
limits given in this paper until it becomes constrained by one of the experi- 
mental points, likely the DIRBE 3.5 /im upper limit. 

Our low and high estimate of the opacity of the universe to TeV and sub- TeV 
7-rays indicate that attenuation should not prevent us from seeing more TeV 
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sources even with existing 7-ray observational facilities. The future projects, 
such as VERITAS, HESS, and MAGIC, have a high potential to increase the 
list of TeV and sub-TeV sources, because they extend their sensitivity range 
for objects similar to the TeV emitter Mrk 501 up to 2; = 0.3. If there are 
more luminous objects capable of producing TeV radiation, this redshift will 
not be a limit. 

The TeV and sub-TeV observations have the potential to further contribute in 
identifying the EBL SED. The HEGRA collaboration has measured recently 
the Mrk 501 spectrum with very high statistical quality (0). For the first 
time the energy range of detected 7-rays is extended above 10 TeV up to 
~ 20 TeV providing a sensitivity in the 40 — 70 /im region of EBL. The 
rapid decline in the Mrk 501 spectrum found at the highest observed energies 
suggests several explanations such as intergalactic EBL absorption due to the 
dust component, intrinsic absorption of 7-rays at the source, or a cutoff in the 
photon production mechanism ([70| ; |69D . Further progress in the analysis of this 
data can be achieved if spectral and temporal characteristics of the emission 
in X-ray and TeV regions are found during multiwavelength campaigns or a 
cutoff feature is identified in the spectrum of Mrk 421 above 10 TeV. In both 
cases, this will reduce the interpretation ambiguity and perhaps the dust EBL 
component expected in this wavelength interval will be detected or strongly 
constrained. 

The valley between 10 — 40 fim may also become testable. The empirically 
based calculations of the EBL spectrum by Malkan & Stecker produce 
a fairly flat EBL SED in this region (see Fig. 10). This feature of their EBL 
model is primarily responsible for a prediction by Stecker (0) of an absorption 
cutoff above ~ 6 TeV for PKS 2155-304 {z = 0.117). It may be that a very 
large exposure on the source may be necessary to deduce spectrum up to this 
energy, but if this cutoff is not observed, it will indicate a SED in this region 
closer to the lines "1-3" shown in Fig. 7 rather than an almost constant EBL 
SED distribution. 

The limits in the wavelength region ~ 1 — 10 /im can be improved if the 
EBL signal in the spectra of AGNs is isolated, or at least limited, by blazar 
emission models. Currently, an upper limit on the spectral curvature of Mrk 
501 and Mrk 421 does not provide a more strict limit than the DIRBE direct 
upper limit at 3.5 /xm. However, if the most likely curvature, 0.17, in their 
spectra is due to EBL absorption then the range of possible SEDs is limited 
by lines marked as "1" and "2" in Fig. 10. These predictions seem to be higher 
than current EBL theoretical models. This could mean either an additional 
contribution to EBL from an unidentified source, or that at least part of the 
blazars' spectral curvature must be intrinsic to the sources. If the latter is 
correct, then confirmation of spectral invariance during episodes of activity 
of blazars may strongly limit flaring mechanisms. Restricting the change of 
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spectral index, ci, due to EBL absorption may impose strict limits on the 
SED in ~ 1 — 10 /im range. Measurements in the GeV range are available 
for Mrk 421 ( [79| ) which show spectral index 1.57. The spectral index for the 
best fit in the 0.26 — 10 TeV range with curvature 0.17 gives value 2.48. The 
0.91 difference again provide a rather loose constraint on EBL since we already 
limited this number to be between 0.12//io and O.GI/Ziq- It is possible to explain 
such a change in the blazar spectral index purely by 7-ray absorption, again 
implying a substantial SED of EBL limited by lines "1" and "2" depending on 
the value of ho. There is no argument, however, to support that such change in 
spectral index could not be intrinsic to the source since these measurements are 
separated by more than two orders of magnitude in energy. The recent report 
of a 5.2a {E > 0.5 GeV) detection of Mrk 501 by EGRET (|0D estimated the 
spectral index as 1.3 ± 0.5 again indicating that the largest allowable opacity 
to TeV photons is possible. It seems that constraints on EBL in the 1 — 10 /im 
wavelength region may come only if the spectrum of a new TeV source will be 
measured, preferably in TeV, sub-TeV, and GeV regions. 

The EBL region 0.1 — 1 fim is not constrained well by current sub- TeV obser- 
vations. This wavelength interval holds crucial information which may finally 
allow one to determine EBL through the whole range ~ 0.1 — 10 /im. VERI- 
TAS, HESS and MAGIC should be able to see farther and detect many sources. 
It is very important that the energy range of these observatories extend to be- 
low 100 GeV providing a possibility to measure in one experiment the "turn 
on" of the effect of sub- TeV 7-ray absorption by EBL. The "knee-like" feature 
has to exist in the spectra of blazars around 100 GeV, independent of proper- 
ties of AGNs. The change of the power law through this knee, if observed for 
several sources, will determine the amount of EBL in this wavelength band. It 
is very important also that due to the high sensitivity of these experiments the 
spectra of the closest AGNs Mrk 421 and Mrk 501 will perhaps be measured 
up to 50 TeV, providing information on the 70 — 140 /xm EBL window. This 
will close the last gap in the EBL SED and meet with direct detections by the 
FIRAS and DIRBE experiments. 
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Fig. 1. Filled squares with arrows indicate EBL upper limits (95%CL) determined 
by various experiments. The wavelength range from 1.25 ^m to 240 fim is from 
the DIRBE experiment (|l8|). Optical and UV bands are by all-sky-photometry ob- 
servations [4l|). The open squares at 140 /im and 240 /im are DIRBE 
detections (shown with la error bars). A tentative detection (|3^ ) is indicated by 
open square marked "T" . The filled circles indicate lower limits (95%CL) obtained 
from galaxy counts. The range 25 — 100 /xm is from IRAS data ([i^). The 6.7 and 15 
fim are from ISO measurements (^). The ground-based measurement at 2.2 ^m is 
reported in (45). The lower limits in the range from 0.36 /im to 0.81 ^m are derived 
from HDF galaxy counts (^6]). The UV lower limit at 0.2;um is from (p7[). The solid 



line is a FIRAS detection of the far- infrared EBL (19). The dashed line indicates 
2.7K CMB radiation. The dotted line is an example of EBL prediction found by 
Primack et al. for the "LCDM Salpeter" model in (|l|). 




Fig. 2. The behavior of the functions f{q), F{q), and q^'^^F{q) which peak at 
q = 0.51, 0.28, and 0.41 respectively. 




Fig. 3. Spectra of Mrk 421 (filled triangles) and Mrk 501 (filled squares) as 
reported in (^). The solid lines indicate the best fits of the curved spectra 
~ exp (^-6\n{E/TeV) - c\n^{E/TeV)^ with parameters {6,c,p^2) for Mrk 421 
(2.50,0.12,0.14), and Mrk 501 (2.26,0.20,0.31). The dotted lines show the best fits 
of the spectra with exponential cutoff ~ exp {—6ln{E/Ec) — E/Ec) and parameters 
{5,Ec,p^2) for Mrk 421 (2.32,7.0 TeV,0.11), and Mrk 501 (1.99,4.9 TeV,0.41). Due 
to variability of the sources the absolute fluxes shown have arbitrary normalization. 



Fig. 4. The chance probabiUty for AGN spectrum to have a certain value of spectral 
index 6. Mrk 421 peak is at = 2.50 with p^2 < for 6 > 2.68, and Mrk 501 

peak is at (5 = 2.26 with p 2 < 10^^ for 6 > 2.44. 




Fig. 5. The chance probabihty for the AGN spectra to have a certain value of 
curvature parameter, c. Mrk 421 peak is at c = 0.12 with p^2 < 10~^ for c > 0.32, 
and Mrk 501 peak is at c = 0.20 with p^2 < 10^^ for c > 0.37. If curvature of 
the spectra is fixed for both AGNs but the spectral indices are allowed to vary, the 
combined chance probability for Mrk 421 and Mrk 501 is given by the dotted line. 
It peaks at c = 0.17 with {p^2 = 0.1), Mrk 421 6 = 2.48, and Mrk 501 6 = 2.28. 
The chance that c > 0.30 is less than 10~^. 
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Fig. 6. Constraints on the choice of pi and p2 parameters from existing EBL upper 
and lower hmits. The po parameter is taken to be 0.025. The dashed hne indicates 
constant attenuation of 10 TeV 7-rays at z = 0.032. The hnes parallel to it would 
correspond to lower or higher constant optical depth. 




Fig. 7. The examples of SED EBL corresponding to maximal and minimal at- 
tenuation of 10 TeV 7-rays consistent with EBL upper and low limits and 
with observations of the Mrk 421 and Mrk 501 TeV spectra. Line 1 de- 
fined by {pi,P2,P3) = (0.025, 0., 0.7) gives minimal optical depth 0.55//io- Line 
3, (0.025, 0., 3.6) generates maximal allowable optical depth 2.88//io- Line 2, 
(0.025,0., 1.6), shows an example of SED EBL which is consistent with recent ten- 
tative EBL detection at 3.5 ^um (36). This EBL distribution would be responsible 
for 1.26/ ho optical depth. 




Fig. 8. The dependence of {co{z), ci{z), C2{z)) functions on redshift for the highest 
allowable EBL absorption effect (marked 1), for the lowest already detected EBL 
density (marked 3), and for the probable EBL density (marked 2, see explanation 
in text). The cq curves are shown as solid lines. The ci and C2 curves are shown as 
dashed and dotted lines respectively. 
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Fig. 9. The change of the observable spectrum of a hypothetical "Mrk 501-like" 
object during its highest activity if it is placed at higher redshifts. Figure "a" shows 
the maximum allowable EBL absorption effect. Figure "b" indicates the opposite 
limit. The sensitivities of Whipple and VERITAS are marked as "W" and "V" 
(see explanation in text). The solid line shows the Mrk 501 spectrum (z = 0.03) 
extrapolated into the 100 — 250 GeV region. The absolute flux of the blazar is 10 
Crab units. The dashed lines correspond to the spectra of such AGN if observed at 
redshift z = 0.06, 0.09, 0.3. 
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Fig. 10. The solid line marked "1" provides an upper limit in the range ~ 1 — 30 
/xm on the SED EBL derived in this work. The range of the EBL SED between 
curves "1" and "2" is preferable if all curvature observed in the spectra of Mrk 
501 and Mrk 421 is attributed to EBL absorption. Upper limits from the Mrk 501 
TeV spectrum by Biller, et al. (28) normalized at the 3.5 /im DIRBE upper limit 
are shown as horizontal bars with arrows. "No cutoff" upper limits from Mrk 421 
spectra are indicated in the range 10 — 40 /xm by solid line segments with arrows (see 
explanation in text). Dotted line is LCDM (Salpeter) model by Primack et al. (IT^), 
and dot-dashed line is LCDM (Scalo) model from the same publication. Dashed line 
is "high" EBL model by Malkan & Stecker (p5|). Filled squares and circles indicate 
upper and low EBL limits respectively (see Fig. 1 for explanation). 



